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Abstract 

The transition density rit and pressure Pt at the inner edge between the hquid core and the solid crust 
of a neutron star are analyzed using the thermodynamical method and the framework of relativistic nuclear 
^ I energy density functionals. Starting from a functional that has been carefully adjusted to experimental binding 

energies of finite nuclei, and varying the density dependence of the corresponding symmetry energy within the 
limits determined by isovector properties of finite nuclei, we estimate the constraints on the core-crust transition 
density and pressure of neutron stars: 0.086 fm"^ < Ut < 0.090 fm"^ and 0.3 MeV fm"^ < Pt < 0.76 MeV fm"^. 
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'=:t 1 Introduction 

O . 

■ Neutron stars are extraordinary astronomical laboratories for the physics of dense neutron-rich nuclear matter 
• • ! [nil]- They consists of several distinct layers: the atmosphere, the surface, the crust and the core. The latter, 
^ ' divided into the outer core and inner core, has a radius of approximately 10 km and contains most of the star's 
k>( , mass. The crust, of w 1 km thickness and containing only a few percent of the total mass, can also be divided into 
; I the outer crust and inner crust. Although less exotic and smaller in size than the core, the crust is nevertheless 
^ \ crucial for the understanding of the physics of neutron stars. It represents the interface between the observable 
surface phenomena and the invisible core. The structure of the crust can be related to some peculiar phenomena, 
such as pulsar glitches, thermal relaxation after matter accretion, quasi periodic oscillations and anisotropic surface 
cooling [31 m [S] . A very important ingredient in the study of the structure and various properties of neutron stars 
is the equation of state (EOS) of neutron-rich nuclear matter 

One of the most important prediction of a given EOS is the location of the inner edge of a neutron star crust. The 
inner crust comprises the region from the density at which neutrons drip-out of nuclei, to the inner edge separating 
the solid crust from the homogeneous liquid core. At the inner edge, in fact, a phase transition occurs from the 
high-density homogeneous matter to the inhomogeneous matter at lower densities. In the transitional region nuclear 
matter exhibits instability against clusterization into a two-phase system: neutron-rich nuclei immersed in dripped 
neutrons (and sometimes protons). As nuclei are arranged in a lattice, they form solid state crust covering the 
star's core, which is considered to be a homogeneous liquid [7]. The uniform matter is nearly pure neutron matter, 
with a proton fraction of a few percent, determined by the condition of /3-equilibrium. The transition density takes 
its critical value ric when the uniform neutron-proton-electron matter (npe) becomes unstable with respect to the 
separation into two coexisting phases (one corresponding to nuclei, the other to a nucleonic sea) 

While the density at which neutrons drip-out of nuclei is rather well determined, the transition density nt at the 
inner edge is much less certain due to our insufficient knowledge of the EOS of neutron-rich nuclear matter. The 
value of nt determines the structure of the inner part of the crust. If sufficiently high, it is possible for non-spherical 
phases, with rod- or plate-like nuclei, to occur before the nuclei dissolve. If nt relatively low, then matter makes a 
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direct transition from spherical nuclei to uniform nucleonic fluid. The extent to which non-spherical phases occur 
will have important consequences for other properties that are determined by the solid crust |8j. 

In general, the determination of the transition density nt itself is a very complicated problem because the 
inner crust may have a very complicated structure. A well established approach is to find the density at which 
the uniform liquid first becomes unstable against small-amplitude density fluctuations, indicating the formation of 
nuclear clusters. This approach includes the dynamical method [SlinilinillllllllllSlIIllIISllIS, the thermodynamical 
method [51 [71 [T71 [TB] , and the random phase approximation (RPA) jTH [50] • 

All theoretical studies have shown that the core-crust transition density and pressure are very sensitive to the 
density dependence of the nuclear matter symmetry energy. The EOS of neutron-rich nuclear matter has been 
constrained by using results from heavy- ion reaction studies [21]. In particular, it has been shown that the Esymip) 
constrained in the same sub-saturation density range as the neutron star crust by the isospin diffusion data in 
heavy-ion collisions at intermediate energies [351 [331 [23] , limits the transition density and pressure to 0.040 fm~'^ 
< Pt < 0.065 fm~'^ and 0.01 MeV/fm^ < -P* < 0.26 MeV/fm^, respectively . These constrained values appear to be 
significantly lower than their fiducial values currently used in the literature. In a very recent study [16j . the core- 
crust transition density and pressure have been systematically analyzed using the dynamical and thermodynamical 
methods with a modified Gogny (MDI) and a set of 51 different Skyrme interactions. Most of these interactions 
predict values for the transition density and pressure that are considerably higher than the intervals cited above. 

In a recent work [25] we have explored a particular class of empirical relativistic nuclear energy density function- 
als, with parameters adjusted to experimental binding energies of a large set of axially deformed nuclei. Starting 
from microscopic nucleon self-energies in nuclear matter, and empirical global properties of the nuclear matter 
equation of state, the coupling parameters of the functional have been determined in a careful comparison of the 
predicted binding energies with data, for a set of 64 axially deformed nuclei in the mass regions A « 150 — 180 
and A « 230 — 250. The isovector channel, in particular, has been carefully adjusted to reproduce available data 
in medium-heavy and heavy nuclei, including neutron-skin thickness and excitation energies of isovector dipole 
resonances. It will be interesting, therefore, to apply this class of relativistic density functionals in a systematic 
investigation of the transition density rit and pressure Pt at the inner edge separating the liquid core from the solid 
crust of neutron stars. In the present study the thermodynamical method will be used. 

In recent years, there has been an increased interest in studies of the relationship between the size of the neutron- 
skin in heavy nuclei, and the symmetry energy at subsaturation densities [TOl [251 [271 [251 [3^ [5011511 15211551 [M1I551 155) . 
Studies have also been reported on the correlation between the size of the neutron-skin and properties of a neutron 
star crust. This was pioneered by Horowitz et. al. [19], who used the Random Phase Approximation based on the 
Relativistic Mean-Field (RMF) framework for nuclear matter and finite nuclei. An almost linear correlation was 
established between the predicted core-crust transition density nt and the size of the neutron-skin. More recently 
such studies have been carried out by by Xu et. al. [Ellll], confirming this linear correlation. 

The article is organized as follows. In Sec. II we review the thermodynamical method used for locating the 
inner edge of a neutron star crust. Sec. Ill contains a brief description of relativistic density functionals that will 
be used to analyze the constraints on the core-crust transition density and pressure of neutron stars. The results 
are presented and discussed in Sec. IV, and Sec. V summarizes the present study. 



2 The Thermodynamical Method 

The core-crust interface corresponds to the phase transition between nuclei and uniform nuclear matter. The 
uniform matter is nearly pure neutron matter, with a proton fraction of just a few percent determined by the 
condition of beta equilibrium. Weak interactions conserve both baryon number and charge |6j , and from the first 
low of thermodynamics, at temperature T — (for details see the Appendix): 

du = —Pdv — fidq, (1) 

where u is the internal energy per baryon, v — 1/n, q = x — Y^, is the charge density, x is the proton fraction in 
baryonic matter, Y^, is the electron density, and in /3-equilibrium ft = Hn — fj,p ^ fi^- The stability of the uniform 
phase requires that u{v,q) is a convex function 37 . This condition leads to the following two constraints for the 
pressure and the chemical potential 

«n-mm>0, (2) 



dvjq \dqj^\dvj^ 
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It is assumed that the total internal energy per baryon u{v, q) can be decomposed into baryon (En) and electron 
(Ee) contributions 

uiv,q) ^ EN{v,q)+Eeiv,q). (4) 
In this work the well know parabolic approximation is used for the baryon energy Eiq{v, q) 

En{v, q) ~ V{v) + Esy,n{v){l - 2xf , (5) 

where higher-order terms in the isospin asymmetry 6 ~ 1 — 2x are neglected. In Eq. ([5]) V{v) denotes the energy 
of symmetric nuclear matter, and Esym{v) is given by 



Esy,n{v) ^ EN{v,q{x = 0)) - EN{v,q{x = 0.5)). 



The electron contribution to the total energy reads 



Ee = -^Ye^J.e 



The condition of /3-equilibrium leads to the relation [381 ES 

'dEr, 



dx 



4:E,ym{n){l - 2x) = ?ic(37r^rie) 



From the relation 



and Eq. (gT]), it follows that 



q — X ~ Ye = X ~ Tle/n, Uf, 



dx 



1 / drie 
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.^AA,(„) V'9AA(„) "Vc'A/i,(„) 8Esy„,{n) nh^c^TT^ 9,Esym{n) fi 
The inequality ^ is equivalent to 

The electron contribution to the pressure Pe is a function of the chemical potential fl — fi^ only 



P. 



J: f4_ 

127r2 {hcY 



and thus the inequality (llip can be written as 

dPb 
dv 



> 0, or n 



dPb 
dn 



> 0. 



(6) 
(7) 

(8) 
(9) 

(10) 
(11) 
(12) 

(13) 



In general the baryonic pressure Pb is a function of both n and x, but for a fixed fi (see also Eq. (H7)) x ~ x{n), so 
that P = P(n,x{n)), and therefore 



m 

dn 



Now, considering that /t = fl{n,x), it follows that 



dPb ^ dPb fdx\ 
dn dx \dnj 



djl 



dp. 
dn 



dn 



dp, 
dx 



dx 



dx 
dn 



dp\ fdfiy ^ 

dn J „ \ dx 



and making use of Eq. (|47)) . one obtains 



dx 
dn 



d^Epf\ ( d'^Ejs! 



Eq. (|26|) now reads 



dPb 
dn 



dndx I \ dx 



2dEN d'^EM /d'^EMY /d'^E 



dn dn2 V dnda 



dx^ 
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The condition of charge neutrality q = requires that x = Yg. This is the case we will consider in the present 
study. After some algebra, it can be shown that the conditions of Eqs. (|17l) and ([TU)) are equivalent to 



Clin) ^n^^+ 2n^ + (1 - 2x)^ 
an'' an 



Ciiin) 



n h 2n — 



An^ 



1 



dn 



E. 



sym 



dE. 



sym 



dn 



>o, 



+ 



> 0. 



(18) 



(19) 



dfij^ 8E 

The second inequality (jl9p is usually valid. The transition density nt is determined from the first inequality ([T| 
For a given EOS, the quantity Ci{n) is plotted as a function of the baryonic density n, and the equation Ci{nt) — 
defines the transition density nt. 



The quantities L and 



sym 



The density dependence of the nuclear matter symmetry energy can be characterized in terms of a few bulk 
parameters by expanding it in Taylor series around the saturation density no 



Esym{n) — Esym{no) + L 



n- no\ Ksym f n - no 



3no 



3no 



0(3)... 



where -Est,™ (no ) = 0-4 is the value of the symmetry energy at saturation, L is the slope parameter 



3no 



\ dn 



and the curvature paraiaeterKsym is the isovector correction to the compression modulus 



K — Qtj^ 



d'^Esynijn) 
dn^ 



(20) 



(21) 



(22) 



The slope parameter L is related to Po , the pressure from the symmetry energy for pure neutron matter at saturation 
density [ID]. The symmetry pressure Po provides the dominant baryonic contribution to the pressure in neutron 
stars at saturation density. It will be interesting to study the relation between the transition density nt , and L and 
Ksym, as well as to examine the correlations between L and Ksym and the neutron-skin thickness. 

The neutron-skin thickness S' of a nucleus is defined as the difference between the root-mean-square radii of 



neutron y^rj) and proton J (r^) distributions 



S=^)-J{rl)=R,,-Rj, 



(23) 



S is sensitive to the density dependence of the nuclear symmetry energy, particularly the slope parameter L 
p6l [3TJ [32] ■ More specifically, the slope parameter L has been found to correlate linearly with the neutron-skin 
thickness of heavy nuclei [TOl [501 [5T| |M] . 



The pressure at the inner edge of a neutron star crust 

The pressure at the inner edge is an important quantity directly related to the crustal fraction of the moment of 
inertia, which can be measured indirectly from observations of pulsars glitches [6] . The total pressure is decomposed 
into baryon and lepton contributions 

P{n,x) = Pb{n,x) + Pe{n,x), 

where 



,dE 



N 



Pb{n, x) = r? 
The baryon pressure P;,, therefore, is given by 

Phin, x) = r? 



EN{n,x) = V{n) + (1 - 2xfEsym{n). 



dV{n) ^ dEsym _ ^^^2 



dn dn 

The electrons are considered as a non-interacting Fermi gas. Their contribution to the total pressure reads 



Pe(n,a;) 



1 



he 



127r2 (nc)3 Un^ 
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(24) 
(25) 



(26) 



(27) 
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For a given symmetry energy Egymin), equation ([8]) determines the equilibrium proton fraction x{n). A simple 
algebra leads to 

x{n) = i - i ([2/3(7 + 1)]'/' - [2/3(7 - 1)]'/') , (28) 



1/2 



1 1 

x(n) — 

where 

/3(n) = 3Tr''n{hc/4Esym{n)j^, l{n) = (^^ + ^ 

3 Relativistic Energy Density Functionals 

The framework of nuclear energy density functionals (NEDF) provides, at present, the most complete microscopic 
approach to the rich variety of structure phenomena in medium-heavy and heavy complex nuclei, including regions 
of the nuclide chart far from the valley of /3-stability. By employing global functionals parameterized by a set of 
w 10 coupling constants, the current generation of EDF-based models has achieved a high level of accuracy in the 
description of ground states and properties of excited states, exotic unstable nuclei, and even nuclear systems at 
the nucleon drip-lines. 

There are important advantages in using relativistic density functionals, i.e. functionals with manifest covari- 
ance. The most obvious is the natural inclusion of the nucleon spin degree of freedom, and the resulting nuclear 
spin-orbit potential which emerges automatically with the empirical strength in a covariant formulation. The con- 
sistent treatment of large isoscalar, Lorentz scalar and vector self-energies provides a unique parametrization of 
time-odd components of the nuclear mean-field, i.e. nucleon currents. Recently |25j we have explored a particular 
class of relativistic nuclear energy density functionals in which only nucleon degrees of freedom are explicitly used 
in the construction of effective interaction terms. Short-distance correlations, as well as intermediate and long- 
range dynamics, are effectively taken into account in the nucleon-density dependence of the strength functionals of 
second-order contact interactions in an effective Lagrangian. 

The basic building blocks are the densities and currents bilinear in the Dirac spinor field ip of the nucleon: 
ipOr^ipj with Ot 6 {IjTi} and T € {1, 7/j,, 75, 757/j, cr^jy}- Here are the isospin Pauli matrices and F generically 
denotes the Dirac matrices. The nuclear ground-state density and energy are determined by the self-consistent 
solution of relativistic linear single-nucleon Kohn-Sham equations. To derive those equations it is useful to construct 
an interaction Lagrangian with four-fermion (contact) interaction terms in the various isospace-space channels: 
isoscalar-scalar (iJip)'^, isoscalar-vector ('tpj^'ilj){'ip'-f^'ip), isovector-scalar {'tpf^p) ■ (V'f-0), isovector-vector {{[jf^^ip) ■ 
{ij;f"ff^tjj) . A general Lagrangian can be written as a power series in the currents ipOrTip and their derivatives, with 
higher-order terms representing in-medium many-body correlations. The Lagrangian considered in Ref. |25j includes 
second-order interaction terms, with many-body correlations (short-distance correlations, as well as intermediate 
and long-range dynamics), encoded in density-dependent coupling functions: 

C ~ iplij ■ d — to)V' 

- ^as(^)('0V')('0'0) - ^"y(^)(V'7''V')(^7MV') - ^o:Tv{n){ipTj^'ip){ipTjf,'ip) 

- ^Ssid.^i^W^i^) - ei,^ ■ ~ ""'V ■ (29) 

In addition to the free-nucleon Lagrangian and the point-coupling interaction terms, when applied to nuclei, the 
model must include the coupling of the protons to the electromagnetic field. The derivative term in Eq. (|29p accounts 
for leading effects of finite-range interactions that are crucial for a quantitative description of nuclear density 
distribution, e.g. nuclear radii. Eq. ([29l) includes only one isovector term, i.e. the isovector-vector interaction 
because, although the isovector strength has a relatively well-defined value, the distribution between the scalar and 
vector channels is not determined by ground-state data. 

The strength and density dependence of the interaction terms of the Lagrangian Eq. ()29|) are parameterized as 
follows: 

as{n) = as + {bs + csx)e-'^'''' , 

av{n) = ay + feye"'^^^, (30) 

where x = n/np, and uq denotes the nucleon density at saturation in symmetric nuclear matter. The set of 10 
parameters has been adjusted in a multistep parameter fit exclusively to the experimental masses of 64 axially 
deformed nuclei in the mass regions A w 150 — 180 and A « 230 — 250. The resulting best-fit functional DD-PCl 
has been further tested in calculations of binding energies, charge radii, deformation parameters, neutron-skin 
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thickness, and excitation energies of giant monopole and dipole resonances. The nuclear matter equation of state, 
corresponding to DD-PCl, is characterized by the following properties at the saturation point: nucleon density 
Psat = 0.152 fm~^, volume energy = —16.06 MeV, symmetry energy 04 = 33 MeV, and the nuclear matter 
compression modulus Knm — 230 MeV. 

4 Results and Discussions 

To adjust the functional DD-PCl, in Ref. [IS] sets of effective interactions with different values of the volume a^, 
surface Os, and symmetry energy in nuclear matter were generated, and the corresponding binding energies of 
deformed nuclei with A w 150 — 180 and A « 230 — 250 were analyzed. The nuclear matter saturation density, the 
Dirac mass, and the compression modulus, were kept fixed: no = 0.152 fm^'' in accordance with values predicted by 
most modern relativistic mean-field models, = 0.58m in the narrow interval of values allowed by the empirical 
energy spacings between spin-orbit partner states in finite nuclei, and Knm = 230 MeV to reproduce experimental 
excitation energies of isoscalar giant monopole resonances in relativistic (Q)RPA calculations. 

Nuclear structure data do not constrain the nuclear matter EOS at high nucleon densities. Therefore, two 
additional points on the E[p) curve in symmetric matter were fixed to the microscopic EoS of Akmal, Pandharipande 
and Ravenhall [41] , based on the Argonne Vig NN potential and the UIX three- nucleon interaction. This EOS has 
extensively been used in studies of high-density nucleon matter and neutron stars. At almost four times nuclear 
matter saturation density, the point n = 0.56 fm~^ with E/A — 34.39 MeV was chosen and, to have an overall 
consistency, one point at low density: n = 0.04 fm~^ with E/A = -6.48 MeV (cf. Table VI of Ref. [IT]). 

The calculated binding energies of finite nuclei are very sensitive to the choice of the nuclear matter volume 
energy coefficient a„. In fact, one of the important results of analysis of deviations between calculated and exper- 
imental masses (mass residuals) of Ref. |25j . is the pronounced isospin and mass dependence of the residuals on 
the nuclear matter volume energy at saturation. To reduce the absolute mass residuals to less than 1 MeV, and to 
contain their mass and isotopic dependence, had to be constrained to a narrow interval of values: —16.04 MeV 

< < —16.08 MeV. Experimental masses do not place very strict constraints on the parameters of the expansion 
of Esym{n) (cf. Eq. ([20]) ). but self-consistent mean- field calculations show that binding energies can restrict the 
values of Egym at nucleon densities somewhat below saturation density, i.e. at n ~ 0.1 fm~'^. Additional informa- 
tion on the symmetry energy can be obtained from data on neutron-skin thickness and excitation energies of giant 
dipole resonances. Recent studies have shown that relativistic effective interactions with volume asymmetry 04 in 
the range 31 MeV < 04 < 35 MeV predict values for neutron-skin thickness that are consistent with data, and 
reproduce experimental excitation energies of isovector giant dipole resonance |42j . Therefore, in the construction 
of the functional DD-PCl in Ref. [25], the volume asymmetry was held fixed at 04 = 33 MeV, and the symmetry 
energy at a density that corresponds to an average nucleon density in finite nuclei: (n) = 0.12 fm~^ was varied. 
The quantity Esymin — 0.12 fm^'^) will be denoted {82)- 

Starting from the relativistic energy density functional DD-PCl, in this work we examine the sensitivity of the 
core-crust transition density nt and pressure Pt of neutron stars, on the density dependence of the corresponding 
symmetry energy of nucleonic matter. In Ref. |25) the value of (5*2) was varied in a rather narrow interval of values 
27.6 MeV < (5*2) < 28.6 MeV, constrained by the empirical values of binding energies and ground-state isovector 
properties of finite nuclei. Fig. [1] displays the corresponding symmetry energy curves Esym as a function of the 
baryon density n. For = —16.06 meV (DD-PCl) the minimum deviation of the theoretical binding energies 
from data is obtained when (5*2) = 27.8 MeV. 

Table [1] and Fig. [2] display the values of the transition density nt (in fm^'^) and transition pressure Pt (in 
Mev-fm"'^), calculated in the thermodynamical model, as functions of < ^2 > for three values of the nuclear 
matter volume energy coefficient a„ . For a given value of the parameter , the values of nt rise with increasing 

< 5*2 >, whereas the opposite is found for the values of Pt. For the considered interval of < 6*2 >, however, the 
changes are small. An increase of 3.5% in < 5*2 > leads to an increase of 1.5% in the value of nt- The transition 
pressure exhibits a somewhat more pronounced dependence (the corresponding decrease is around 16 — 20%). Both 
nt and Pt display a negligible dependence on a^, even though a^, = —16.02 MeV and a„ = —16.14 MeV lie outside 
the interval of values for which the absolute deviations between calculated and experimental masses are smaller 
than 1 MeV. 

In Fig. |3|we plot the transition pressure Pt as a function of the transition density nt for the three sets of nuclear 
matter EOS and symmetry energy described above, in comparison with results of recent calculations performed 
using an isospin and momentum-dependent modified Gogny effective interaction (MDI) [16l [43]. The different 
values of the parameter x in the MDI model correspond to various choices of the density dependence of the nuclear 
symmetry energy. In Refs. |241 144] it has been shown that only — 1 < a; < leads to a density dependence of 
the symmetry energy in the sub-saturation density region that is consistent with isospin diffusion data and the 
empirical value of the neutron-skin thickness in ^"^Pb. In addition to the MDI EOS, in Fig. [3] we also show the 
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result obtained by Akmal et al. f4P with the A18 + Sv + UIX* interaction (ARP), and the value obtained in 
the recent Dirac-Brueckner-Hartree-Fock (DBHF) calculation ^32j with the Bonn B One-Boson-Exchange (OBE) 
potential (DBHF-hBonn B) 

A distinctive feature of the present analysis is the narrow interval of allowed values {nt,Pt) that results from 
the rather stringent constraints on the parameter < 5*2 >. The effect of varying the volume energy at saturation 
is almost negligible. The present results for rit and Pf lie in the region constrained by the measure of the current 
uncertainty in the density dependence of the symmetry energy 6 , and are found very close to the result of Akmal 
et al. [H]. We note that all the results shown in Fig. [3] are obtained using the parabolic approximation for the 
EOS of isospin-asymmetric nuclear matter. The transition density and pressure have also been estimated using 
the full equation of state and employing both the dynamical and thermodynamical methods |151 116j . 

As explained above, the rather narrow interval of < 5*2 >, for this type of nuclear energy density functionals, 
has been constrained by the empirical values of binding energies and ground-state isovector properties of finite 
nuclei. The symmetry energy at saturation density, 04 = 33 MeV, was fixed in Ref. to obtain the best 
results for the neutron-skin thickness in Sn isotopes and ^°®Pb, and for the excitation energies of isovector dipole 
resonances. However, because of large experimental uncertainties, especially for the neutron-skin thickness, good 
agreement with data can also be obtained for other values of 04. This is shown in Fig. |4l where we plot the 
predictions for the differences between neutron and proton rms radii of Sn and Pb isotopes, in comparison with 
available data [151 SZl HH] , for different choices of the symmetry energy at saturation density. The self-consistent 
mean- field calculations have been performed using the relativistic Hartree-Bogoliubov (RHB) model [35^, with 
pairing correlations described by the pairing part of the finite-range Gogny interaction. The isoscalar channel of 
the particle-hole interaction corresponds to the DD-PCl functional, and in the isovector channel < 6*2 > is kept 
fixed at 27.8 MeV (DD-PCl), whereas 04 is varied in the interval between 30 MeV and 35 MeV. The corresponding 
symmetry energy as a function of the nucleon density is shown in Fig. O 

We notice, therefore, that by keeping < 5*2 > constant and varying a4 in the interval between 30 MeV and 35 
MeV, the density dependence of the symmetry energy can be modified in a controlled way, i.e. the corresponding 
energy density functionals still reproduce ground-state properties of finite nuclei in fair agreement with data. In 
Table [2] and Fig. [6] we display the corresponding values of the transition density nj (in fm^'^) and transition 
pressure Pt (in Mev-fm"'^) as functions of 04 for three values of the nuclear matter volume energy coefficient 
The transition pressure Pt as a function of the transition density for the three sets of nuclear matter EOS 
and symmetry energy is plotted in Fig. [T] Not surprising, considering the symmetry energy curves of Fig. [5l the 
constraints on nt and Pt have been relaxed in this case, and the allowed values span a much larger interval of values 
compared to the restricted variation of < 6*2 > shown in Figs. [l]and[2j 

To be able to compare the present results for the transition density and transition pressure with recent studies 
|16) . in Fig. [5] we plot the calculated values oirit and Pt as functions of the slope parameter of the symmetry energy 
(cf. Eq. (HIJ), for the two sets of effective interactions described above, nt is a monotonously decreasing, and Pt 
monotonously increasing function of L. In the small interval of L values determined by the variation of < 6*2 > 
between 27.6 MeV and 28.6 MeV, both nt and Pt display a linear dependence on L. In the much larger interval 
determined by the variation of 04 from 30 MeV to 35 MeV, a weak parabolic dependence of nt and Pt is found. We 
note that transport model studies of the isospin diffusion data in heavy-ion reactions have constrained the slope 
parameter L to the values 88 ±25 MeV [15]. Considering that we can also, most probably, exclude the value 04 = 30 
MeV for the asymmetry at saturation density (cf. Figs. |4] and [6]), because it implies an unrealistically small value 
of < 0.1 fm for the neutron-skin thickness of ^"^Pb, the present analysis places the following constraints on the 
core-crust transition density and pressure of neutron stars: 0.086 fm~'^ < nt < 0.090 fm~^ and 0.3 MeV fm~'^ < 
Pt < 0.76 MeV fm~^ 

Finally, in Fig.|S]we compare the present prediction for the range of values of the transition density nt with the 
results of Horowitz and Piekarewicz who, in Ref. [19) . also used the framework of relativistic mean- field effective 
interactions to study the relationship between the neutron-skin thickness of a heavy nucleus and the properties 
of neutron star crusts. Starting from the NL3 meson-exchange effective interaction |50j . the density dependence 
of the symmetry energy was varied by adding nonlinear couplings between the isoscalar and isovector mesons 
to the original interaction. The variation was carried out in such a way to enhance the changes in the neutron 
density and neutron-skin thickness, while keeping small the corresponding changes in the binding energy and proton 
density distribution. For the solid crust of a neutron star, the effective RMF interactions were used in a simple 
RPA calculation of the transition density below which uniform neutron-rich matter becomes unstable against small 
amplitude density fluctuations. The resulting transition densities are plotted in Fig.|9]as a function of the predicted 
difference between neutron and proton rms radii in ^"^Pb. This inverse correlation was parameterized [19j 

rif w 0.16 - 0.39(i?„ - iip), (31) 

with the skin thickness expressed in fm. In the present analysis, using a different type of relativistic effective 
interactions and varying the density dependence of the symmetry energy by explicitly modifying < ^2 > or 04, we 
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find a much weaker dependence nt on the neutron-skin thickness of ^"^^Pb. 

5 Summary 

The framework of relativistic nuclear energy functionals has been employed to analyze and constrain the transition 
density rit and pressure Pt at the inner edge between the liquid core and the solid crust of a neutron star, using the 
thermodynamical method. Starting from a class of energy density functionals carefully adjusted to experimental 
masses of finite nuclei, we have examined the sensitivity of the core-crust transition density rit and pressure Pt 
on the density dependence of corresponding symmetry energy of nucleonic matter. The limits of variation of the 
density dependence of the symmetry energy are determined by isovector properties of finite nuclei: the thickness 
of the neutron-skin and the excitation energies of isovector giant dipole resonances. Instead of an unrestricted 
variation of the parameters of the Taylor expansion of the symmetry energy around the saturation density of 
nuclear matter, that is the slope parameter and the isovector correction to the compression modulus, we modify 
the density dependence by varying the value of the nuclear symmetry energy at a point somewhat below the 
saturation density {S2) (the symmetry energy at n = 0.12fm~'^), and at the saturation density a4 (the symmetry 
energy at n = 0.152fm~"^, the saturation density for this class of relativistic density functionals). In the former 
case, for a given value of the volume energy coefficient a^, {S2) has been varied in a rather narrow interval of 
values 27.6 MeV < {S2) < 28.6 MeV determined by a fit to the experimental binding energies. We have found 
that an increase of 3.5% in < 5*2 > leads to an increase of 1.5% in the value of nt while Pt exhibits a somewhat 
more pronounced dependence (the corresponding decrease is around 16 — 20%). Both nt and Pt display a negligible 
dependence on ay. The variation of the parameter 04 has been in the range of values: 30 MeV < 04 < 35 MeV, 
allowed by the empirical thickness of the neutron-skin and excitation energies of isovector dipole resonances, for a 
fixed value of < 6*2 >. Again, there is virtually no dependence on ay, but now both nt and Pt span much wider 
intervals. We have also examined the dependence of nt and Pt on the slope parameter of the symmetry energy L, 
for the two sets of effective interactions described above. For the empirical range of the slope parameter 88 ± 25 
MeV, and comparing the calculated values of the neutron-skin thickness with available data for Sn isotopes and 
^"^Pb, we have deduced the following constraints on the core-crust transition density and pressure of neutron stars: 
0.086 fm"^ <nt< 0.090 fm"^ and 0.3 MeV fm"^ < Pt < 0.76 MeV fm"^ 

The present study will be extended to include also the dynamical and the random phase approximation ap- 
proaches to analyze the nuclear constraints on the core-crust transition density and pressure of neutron stars. 
Interesting issues for future study will also be finite temperature and neutrino trapping effects. Work along these 
lines is in progress. 
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Appendix 

Proof of the equality: du = —Pdv — p,dq 

From the first low of thermodynamic, at temperature T = 



where U is the total energy of = Ni particles in a volume V. Here i — p,n, e, and the number of baryons is 
A^fa = A'ji -I- Np. Dividing Eq. ([5^ by the baryon number Nf,, one obtains 




(32) 




(33) 



with 




(34) 
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and 



Y — ^ —Y 



A simple algebra leads to the following relation 

tJ-idYi = -{^J,n - A*p)dYp + ^edYe = -fldYp + fidYe 

i 

= -fi{dYp - dYe) = -fidq, 

where we have used the equalities 

= IJ-n - tJ'p = IJ-e, and q = Yp-Ye. 

Eq. therefore takes the form 

du = —Pdv — /idq. 

Convexity of the function u{v,q) 

Let us consider the function u{v,q) and the determinant 



(35) 

(36) 

(37) 
(38) 



D = 



dv^ dvdq 
dqdv dq'-^ 



The differential du(u, q) reads 



From Eqs. ^ and (gO]) it follows that 



du 



dv 



du{u, g) = — dv + — dq 



du 



du 
dv 



' du 

M = - 1 IT- 



The determinant Eq. (H^ takes the form 



D = 



The requirement of convexity for the function u{v,q) leads to the following two sets of inequalities 
Case 1: 

D>0, - ( ^ I > 0. 



Case 2: 



The second set can written in the form 



D>0, 



dv 



f-\ >0. 



9qJ ^ 



dvjAdqJ^ \dqj^\dv 



>o. 



> 0, 



(39) 
(40) 
(41) 

(42) 

(43) 
(44) 

(45) 
(46) 



dfi 



dv 



dfl — [ — ] dv + [ — ] dq, 



dfi 



Dividing the inequality (|45p by the positive quantity ~ , and considering that the differential of the function 

jji{v,q) is given by 

(47) 
(48) 

> 0. (49) 



and therefore for a fixed {1 



one obtains the following inequality 



djl 



dfi\ f dq\ ^ 



9v ) q p 



dP\ fdP\ f dq 



dv 



dq J \dv 
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Table 1: The values of the transition density rit (in fm^'^) and transition pressure Pt (in Mev-fm~^), calculated 
in the thermodynamical model, as functions of < 52 >, for 04 = 33 MeV and three values of the nuclear matter 
volume energy coefficient a„. 
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0.0868 


0.598 
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Tabic 2: The vahics of the transition density rit (in fm^"^) and transition pressm'e Pi (in Mev-fm^"'), calculated 
in the thermodynamical model, as functions of 04, for < S2 >= 27.8 MeV and three values of the nuclear matter 
volume energy coefficient a„. 
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Figure 1: The symmetry energy E^ym as a function of the nucleon density n for various values of the parameter 
< 5*2 >. The symmetry energy at saturation is 04 = 33 MeV. 




Figure 2: The transition density rit (a), and the transition pressure Pt (b), as functions of < 5*2 > for three values 
of the nuclear matter volume energy coefficient a^. 
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Figure 3: The transition pressure Pt as a function of the transition density n*. For a fixed value of the symmetry 
energy at saturation = 33 MeV, and three values of the nuclear matter volume energy coefficient a„, the 
parameter < 5*2 > is varied in the interval between 27.6 MeV and 28.6 MeV. The resulting constraints of Pt and 
nt are plotted in comparison with results obtained using different models (for details see Ref. [43]). The lines 
Pt = 0.2 MeV fm~'^ and Pt — 0.65 MeV fm~^ correspond to the measure of the current uncertainty in the density 
dependence of the symmetry energy [5] . 
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Figure 4: RHB predictions for the differences between neutron and proton rms radii of (a) Sn and (b) Pb isotopes, 
in comparison with available data [461 1471 148) , for different choices of the symmetry energy at saturation density 
a4 (in MeV). 
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Figure 5: The symmetry energy Esym as a function of the nucleon density n for various values of the symmetry 
energy at saturation a^. The parameter < 5*2 > (see text) is kept constant at 27.8 MeV. 
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Figure 6: The transition density rit (a), and the transition pressure Pt (b), as functions of the symmetry energy at 
saturation density 04, for three values of the nuclear matter volume energy coefficient a^. 
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Figure 7: Same as described in the caption to Fig. |3] but for fixed < 5*2 >= 27.8 MeV, and the symmetry energy 
at saturation in the interval 30 MeV < a4 < 35 MeV. 




Figure 8: The transition density nt (a), and the transition pressure Pt (b), as functions of the symmetry energy at 
saturation density 04, for three values of the nuclear matter volume energy coefficient ay 



16 




0.06 



0.12 



0.16 



0.20 

R -R (fm) 



0.24 



Figure 9: The transition density rit as function of the neutron-skin thickness i?„ — Rp of ^°*Pb. The values of rit 
calculated using the thermodynamical model in the present work (solid) , are compared with those of Ref. [TH] (see 
text for description). 
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